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We present a model of the compact star in Cassiopea A, which accommodates its unusually fast cooling behavior. This feature is inter- 
preted as due to an enhancement in the neutrino emission triggered by a transition from fully gapped two- flavor color superconducting 
phase to a crystalline or an alternative gapless color superconducting phase. By fine-tuning a single parameter - the temperature of 
this transition - a specific cooling scenario can be selected that fits the Cas A data. Such scenario requires as a prerequisite a massive 
M ~ 2M Q star and is, therefore, distinctive from models invoking canonical 1.4 M Q mass star with nucleonic pairing only. 
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1. Introduction 

Recently, an unprecedented fast cooling of the compact star in 
Cassiopeia A (Cas A) - the youngest known supernova remnant 
in Milky Way - was i nferred from the Chandra observations ove r 
a period of 10 years (iHeinke & Holl2010t IShternin et ai]|2011l) . 
The age of the object 330 yr is accurately deduced from the as- 
sociation of the remnant with the historical supernova SN 1680. 
This age estimate is confirmed through kinematical estimates. 
The fits to the surface temperature of the star suggest that the 
thermal soft x-ray spectrum originates from a non-magnetized 
carbon atmosphere at temperature ~ 2. x 10 6 K and emitting 
region of 8-17 km. 

Hadronic models of compact stars have been invoked to fit 
the data by assuming a canonical 1.4M mass neutron star. The 
models proposed include those which attribute the Cas A cooling 
to a star that (a) undergoes a transition to superfluid state in the 
baryonic core, which is accompanied by an en hancement of the 
cooling due to Coop er pair-breaking processes dPage et al.l201 it 
IShternin et al.ll20 1 lb ; (b) experiences slower than commonly ac- 
cepted thenrml^ekxationhi the core due to low thermal conduc- 
tivity dBlaschke et al.ll20T2h : (c) undergoes rotationally induced 
compres sion which triggers a r apid Urea cooling in the center of 
the star dNegreiros et al.ll2013l) . 

Recent observations of PSR J1614-223010+1 1 lend firm 
evidence that equilibrium seq uences of neutron sta rs contain 
massive ~ 2M Q members (iDemorest et al.l l2010h . Models 
which predict quark matter in such massive compact stars 
in col or-superconducting s tate have been st udied exten- 
sively (lAlford et al.l l2005at iKlahn et al l 120071: llppolito et alJ 
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important, because, otherwise, the unpaired quarks neutrino- 
cool a star rapidly to temperatures well below of t hose observed 
~" lAlford et al.ll2005bl: lAnglani et al.l 
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The true ground state of QCD at densities relevant for com- 
pact stars is not known. At relevant densities the strong coupling 
as ~ 1, i.e., perturbation theory is not useful. More importantly, 
because of /^-equilibrium and strange quark mass, cross-flavor 
light, i.e., up- and down-quark pairing occurs among the quarks 



living on different Fermi surfaces. As a consequence, there are 
modifications to the standard Bardeen-Cooper-Schrieffer (BCS) 
picture of pairing. The non-BCS pairing may lead to one of 
the possible gapless two-flavor p hases ( Shovkovv & Huangl 
12003b iMuther & Sedrakianl 120031: IShovkovv et al.l l2005h or. 
for example. 
(lAlford et al.1 
2005 
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Huang & Sedrakianll2010t Unglani et al.ll2013l) . The crvstalline- 
color-superconductivity comes in several variants, which 
differ by the way the translational symmetry is broken by the 
condensate of Cooper pairs carrying finite momentum. Below 
we shall use results obtained for the so-called Fulde-Ferrell 
(hereafter FF) phase, which is simple to model, but is general 
enough to preserve key feature of the crystalline phases, 
namely the existence of gapless modes on the Fermi s u rfaces 
of quarks (lAlford et al. 1200 it ISedrakian & Rischkel 12009b 
iHuang & Sedrakianll20ldl) . 

Proto-neutron stars emerge from supernova explosions with 
temperatures T ~ 50 MeV and nearly isospin symmetric matter 
in the core. After the initial fast cooling to T ~ 0. 1 MeV the core 
reaches its asymptotic state which is highly isospin asymmetric. 
The initial fast temperature drop is followed by much slower 
cooling to temperatures of order 10 keV over the following 
10 4 — 10 5 years. Initially, high-temperature and low-isospin quark 
matter core is in the perfect two-flavor color-superconducting 
(hereafter 2SC) state where the Fermi surfaces of the up- and 
down-quarks are completely gapped. As the star cools and the 
/3-equilibrium is established, the 2SC phase becomes unsta- 
ble, because the asymptotic ground state of quark matter corre- 
sponds to a class of gapless and/or crystalline two-flavor color- 
superconducting phases mentioned above. This is illustrated in 
the top panel of Fig. [1] where we show the phase diagram of 
the FF and 2SC phases i n the temperature and flavor asymme- 
try plane [adapted from ISedrakian & Rischke I (120091) 1. During 
the early evolution of the star, the matter traverses from the up- 
per left corner of the phase diagram (i.e. from unpaired, almost 
isospin symmetric matter) to the right lower corner, which cor- 
responds to the cold FF phase and large isospin asymmetry. The 
phase diagram of the same phases in the density-temperature 
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Fig. 1. Upper panel. The phase diagram of color- 
superconducting quark matter in the temperature (T) and 
isospin imbalance a = (rid - n u )/(n u + plane, where rid and 
n„ are the number densities of d and u flavors of quarks. The 
arrows show schematically the path in the phase diagram during 
the cooling of a star. Lower panel. Phase diagram of the same 
phases, but in the temperature and baryon density (ng) plane for 
matter in /^-equilibrium with electrons. 



plane and under /j-e quilibrium is shown in the lower panel of 
Fig. [J [adapted from lHuang & SedrakiarJ(l2010l) l. 

In this work we conjecture that the rapid cooling of Cas A 
can be understood as a phase transition from the perfect 2SC 
phase to a crystalline/gapless color-superconducting state. Thus, 
somewhat counter-intuitively, the rapid cooling is the result of a 
phase transition within the phase diagram of dense QCD from 
one to another color-superconducting phase. This transition is 
characterized by a temperature T* and the transient "widths" w, 
characterizing the finite time-scale of the transition. 

In general T* will depend on a number of factors, such as 
the nucleation history of the superconducting phase, geometry 
of the superconducting regions, dynamics and order of the phase 
transition, etc. Furthermore, because of the density dependence 
of the critical temperature T* in the stellar interior, this phase 
transition will not occur globally in the quark core, but gradu- 



ally within some shells. Therefore, the values of T* and w will 
be determined by the local conditions in shells. We expect that 
the transition is not instantaneous: it takes place via a propa- 
gation of the phase-separation fronts rather than a coherent in- 
stantaneous transition. The equilibrium model calculations give 
T c ~ 10 MeV in a particular (Nambu-Jona-Lasinio Lagrangian 
based) model, see Fig. Q] But the absolute value of the gap is 
strongly model dependent and it may differ from the equilibrium 
critical temperature of phase transition for dynamical reasons. 
(Supercooled state in electronic superconductor are observed, 
where the system remains in the initial unpaired state below the 
equilibrium critical transition temperature T c , because of the fi- 
nite time-scales of the phase transition process). These consid- 
erations motivate our treatment of T* as a free parameter which 
will be used to fit our cooling model to the Cas A data. These 
fits are not sensitive to the precise value of transient width w 
characterizing the duration of the transition. 

A preliminary repo rt of the present study was given else- 
where (ISedrakianl2013h . 



2. Modeling cooling processes 

Our cooling model is described in iHess & Sedrakianl d20TTl) . 
where the details of the microscopic input and the method of 
solution can be found. Here we describe the extension of this 
model, which allows us to accommodate the phase-transition be- 
tween the phases in a phenomenological manner. 

Quark matter consisting of two light u and d flavors of 
quarks cools predominantly via the beta decay (Urea) reactions 
d —* u + e + v, u + e — > d + v, where e stands for electron, v 
and v represent electron neutrino and antineutrino. For unpaired 
quarks and to leading order in the strong coupling constant as 
th e emissivi t y of t he process per quark color was calculated 
bv llwamotold 19801) 

£p = cos2 0pdP„p e a s T 6 , (1) 

where G is the weak coupling constant, 8 is the Cabibbo angle, 
Pd, Pu, and p e are the Fermi momenta of down quarks, up quarks, 
and electrons. 

Quark pairing modifies the temperature dependence of the 
process (Q]i. In the BCS-type superconductors, the process is 
suppressed almost linearly for T - T c and exponentially for 
T <k T c , where T c is the critical temperature. In gapless su- 
perconductors the emergence of the new scale dp = (pd - M«)/2> 
where p u j are the chemical potentials of light quarks, leads to 
two qualitatively different possibilities, which are distinguished 
by the value of the dimensionless parameter t — Ap/du , where 
Ao is the gap in the limit dp - ( Jai kumar et al.ll2006b . When 
£ > 1, the entire Fermi surface is gapped. Consequently, when 
the thermal smearing of the Fermi surface ~ T is much smaller 
than the gap, no excitation can be created out of the Fermi 
sphere. When f < 1, particles can be excited from the gapless 
regions of the Fermi sphere. For the FF phases the shift in the 
chemical potential is replaced by the anti-symmetric in flavor 
part of the single particle spectrum of up- and down-quarks, e u /d, 
i.e., dp — > [ed(Q) - c«(Q)]/2, where Q is the total momentum of 
a Cooper pair. In the case of single plane wave crystalline phase 
the spectrum depends on the magnitude of the vector Q only. 
The asymptotic behavior of the f parameter discussed above, 
permits us to use an interpolation formula, which covers the 
asymptotic limits, and does not require detailed modeling of the 
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Fig. 2. Left panel. Dependence of the surface temperature on time for the values of the transition temperature T* (in MeV) indicated 
in the plot and for fixed width w = 0.2 MeV. The underlying model is the 1.91 M Q compact star model of iKnippel & Sedrakianl 
d2009t) . Right panel. Dependence of the surface temperature (in units of 10 6 K) on time for the same model. The points with error 
bars show the Cas A data, the solid lines are fits to these data by variation of T* (again for fixed width w = 0.2 MeV). 



physics of the underlying phase dJaikumar et al.ll2006h 

T<T C ) = 2/(0*. f(0 = 



exp[(£ ■ 



1] + 1 



(2) 



where the factor two is the number of colors. 

Because the 2SC pairing pattern breaks the SU(3) color sym- 
metry, one of the quark colors (say, the blue color) is not in- 
volved in that pairing. The strength and the flav or content of 
pairing among blue quark s is model de pendent (lAlford et al.l 
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2006; Aguilera 2007). It is likely that cross-flavor (blue-up and 
blue-down) pairing is suppressed by large mismatch between 
the Fermi surfaces of the up and down quarks. Pairing is ex- 
pected in the color 6s and flavor 3$ channel, which is same- 
flavor and same-color pairing, and therefo re is not affected by 
the flavor asymmetry (lAlford et al.l 120031) . However, the pair- 
ing in this channel vanishes in the chiral limit, when the mass 
of blue quarks vanishes. Pairing due to alternative mechanisms, 
outside those studied so far is not excluded. Quite generally, sin- 
gle flavor and/or color pairings have gaps in the range of 10-100 
keV, which are much smaller than the gaps emerging from 2SC- 
type pairing. Since there are no definitive theoretical prediction 
about the pairing in the blue sector we will assume that the neu- 
trino emissivity of blue quarks is suppressed exponentially by 
the pairing gap, as is the case for nucleonic matter. Therefore, 
we shall simply parameterize the modifications of the Urea pro- 
cess on blue quarks as in the case of baryonic matter, i.e., 



^(r<r c ) = 6*ex P 



Ab 
' T 



(3) 



where At, is the pairing gap for the blue color quarks. Note 
that e b p = €p 12 due to the different number of colors involved. 



Several studies used similar p arameterizations dGrigorian et al.l 
120051 : IHess & Sedrakianll26TlT) . Furthermore, we will choose to 
work with a pairing gap A/, = 100 keV, which has the conse- 
quence that the blue quarks do not affect the cooling dynamics 
and the fits to Cas A data, i.e., they play only a passive role. The 
effect of s maller values of this gap on the cooling is discussed 
elsewhere (IHess & Sedrakianll201 ll) . 

In addition to neutrino emissivities we need to model the 
specific heat of the quark phases. For the inhomogeneous 
red -green condensate the critical temperature changes with 
6n (IHess & Sedrakianll201 lb 



T c (0 - TVo 



(4) 



where fi = fa + M „)/2, A = A(£ = 0), T c0 = T c (( = 0), 
and 5(0 = {rid - n„)/(rid + n u ). The fully gapped and gapless 
regimes behave differently; in the presence of gapless modes, 
the specific heat has a linear dependence on the temperature. A 
phenomenological way to model this behavior is given by the 
relation 



c r s g (C,T<T c ) = f(Oc r , 



(5) 



where c' g is the specific heat of red-green unpaired quarks, taken 
as that for noninteracting quarks and cj? is the specific heat of 
pair-correlated quarks. To compute the specific heat of the elec- 
trons, we assume a noninteracting, ultrarelativistic gas. 

Finally, we need the relation between the temperature of 
the isothermal core and the crust, which can be written as 
dPotekhin et alJI 1997b 

T s e = (or 8 /*J /4 (6) 

where T s ^ is the surface temperature in units of 10 6 K, and T% 
is the isothermal core temperature in units of 10 8 K, a and /3 are 
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constants which depend on the composition of the star's atmo- 
sphere, g s is the surface gravity in units of 10 14 cm s~ 2 . We will 
adopt the values /3 = 0.55, which is intermediate between the 
purely-iron (/J = 0.5) and fully accreted envelope (J3 = 0.61) and 
a = 181 appropriate for an accreted envelop. The fit results are 
sensitive to the first parameter, but not the second. 

3. Simulation results 

The key idea of our model is that the transition occurs from the 
fully gap (perfect) 2SC phase to some type of crystalline/gapless 
phase. A particular realization of this scenario would be a tran- 
sition to single plane wave crystalline phase (the FF phase). We 
model this transition phenomenologically by adopting tempera- 
ture depe ndent (and, therefore, time -dependent) parameter f (T). 
[Note that lHess & Sedrakianl(l2011l) assumed constant f and thus 
phase transitions of the type discussed here were excluded in 
their study from the outset]. The temperature/time dependence 
is modelled such that for T > T* one obtains £ > 1, whereas for 
T < T", £ < 1. The £(T) function is modelled by the simply law 

aT) = St-^g(T), (7) 

where £ is the initial value, A£ is the constant change in this 
function, and the function g(T) affects the transition from the 
initial value £ to the asymptotic final value = - S(. The 
transition is conveniently modelled by the following function 

g(T) = . (8) 

exp(I^)+l 

which allows us to control the temperature of transition, fixed by 
7**, and the smoothness of the transition, fixed by the width w. 

Fig. |2 left panel, shows the cooling evolution of M — 
1.91M© model of a compact star for several values of T* and 
fixed w = 0.2 MeV. It is seen that for a specific value of T* the 
cooling curve passes through the location of Cas A in the log 
7*-log t plane. The right panel of Fig. |2]shows the observed evo- 
lution of Cas A for the same value of w and two (fine-tuned) 
values of T*; the best fit to data is shown by solid line; the 
dashed line demonstrates the sensitivity of the fit to variations 
of the parameter T*. In each case the critical temperature of 
the same-color pairing of the blue quarks satisfies the condition 
T c b = A;,/ 1.76 > T*, which implies that blue quarks have no 
effect on the dynamics of the transition process. The numerical 
value A/, = 0.1 MeV was used in the simulations. In case of low 
critical temperatures of blue quark pairing a compact star will 
rapidly neutrino-cool on time-scales t < 300 yr below the tem- 
perature of Cas A, which obviously will preclude any fit to the 
data. 

4. Conclusions 

We have demonstrated that the rapid cooling of Cas A can be 
understood as a phase transition from the perfect two-flavor 
color-superconducting phase to a gapless/crystalline phase. 
Counterintuitively, the phase transition occurs within the phase 
diagram of dense QCD and reflects the ordering of various su- 
perconducting phases in the temperature, density and isospin 
spaces. This idea is supported qualitatively by the structure of 
the phase diagram of the dense QCD at densities and tempera- 
tures relevant to cooling evolution of compact stars. 

We have shown that the data can be fitted by adjusting a sin- 
gle parameter - the phase transition temperature - at a fixed value 



of transient time-scale, which encodes the dynamics of the tran- 
sition. A further important assumption is that there is no interfer- 
ence from the blue quarks, so that the temperature of the model 
is above the Cas A temperature for t < 300 yr, which is an ob- 
vious requirement for any successful fit to data in the absence of 
internal heating. 

The present model implies that only massive compact stars 
undergo a rapid temperature drop observed in Cas A, for only 
massive members of compact star sequences contain (color- 
superconducting) quark matter in their interiors. This feature 
sets our model apart from the alternatives which invoke nu- 
cleonic pairing and pair-breaking processes dPage et al.l 120 1 U 
I Shternin et al.ll2.01 lb or slow thermal relaxation dBlaschke et al.l 
120121) . i.e., physics of canonical mass neutron stars. Rotation 
is unimportant for ou r modeling as opposed to the model of 
iNegreiros et al.l (120131) . Clearly, the various models can be dis- 
criminated through their distinctive features in case further ob- 
servational information will become available on the properties 
of Cas A. Evidently, it would be premature to conclude that the 
Cas A behavior is a direct evidence for superfluidity or even mea- 
surement ofT c for any particular type of a superfluid in neutron 
star interiors. 
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